Conventional character of the BCS-BEC cross-over in ultra-cold gases of '^"K 
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We use the standard fermionic and boson-fermion Hamiltonians to study the BCS-BEC cross-over near the 
202 G resonance in a two-component mixture of fermionic ''"K atoms employed in the experiment of C.A. Regal 
et al. [Phys. Rev. Lett. 92, 040403 (2004)]. Our mean-field analysis of many-body equilibrium quantities shows 
virtually no differences between the predictions of the two approaches, provided they are both implemented 
in a manner that properly includes the effect of the highest excited bound state of the background scattering 
potential, rather than just the magnetic-field dependence of the scattering length. Consequently, we rule out 
the macroscopic occupation of the molecular field as a mechanism behind the fermionic pair condensation and 
show that the BCS-BEC cross-over in ultra-cold """K gases can be analysed and understood on the same basis as 
in the conventional systems of solid state physics. 

PACS numbers: 03.75.Ss, 03.75.Nt, 34.50.-.S 



I. INTRODUCTION 

Dilute vapours of fermionic atoms have recently attracted 
considerable attention as systems for studying the cross-over 
between Bardeen-Cooper-Schrieffer (BCS) pairing and Bose- 
Einstein condensation (BEC) of self-bound pairs of particles 
ill, i2| |3|, |4j, |5|1 . The experimental techniques underlying these 
studies all take advantage of the unique possibility of con- 
trolling the inter-atomic interactions via magnetically tunable 
Feshbach resonances. Observations of both molecular BEC 
ll^QI^ and condensation of unbound fermionic pairs llgllT^ 
in gases containing incoherent mixtures of two different spin 
components have recently been reported. The properties of 
both phases have been studied in several recent experiments 
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The principal question of whether the BCS-BEC cross-over 
in ultra-cold atomic gases can be analysed and understood 
on the same basis as in the traditional systems of solid state 
physics d H H H is, however, a matter of continuing 
controversy. This controversy dates back to the first predic- 
tions of a superfluid phase transition in dilute Fermi gases us- 
ing Feshbach resonances iUtI ITsIi . suggesting a new type of 
pair condensation based on a macroscopic occupation of the 
Feshbach resonance level. The concept underlying this idea 
of resonance superfluidity iUtI ITsIl was motivated by the use 
of an effective boson-fermion Hamiltonian JTgl Eoll . which 
treats pairs of atoms in the resonance state configuration in 
terms of single structureless Bose particles, often referred to 
as "molecules". The studies of Refs. 1I21I l22il have clearly 
shown, however, that the experimentally relevant highest ex- 
cited vibrational multi-channel molecular bound state (the 
Feshbach molecule) is, in general, significantly different from 
the Feshbach resonance level (belonging to the closed scatter- 
ing channel) and requires an explicit description in terms of a 
composite two-body system. 

In this paper we demonstrate that the thermal equilibrium 
physical quantities, relevant to recent experiments on the 
BCS-BEC cross-over in cold gases using broad Feshbach res- 
onances of '**'K and ^Li, can all be described by the usual 
fermionic Hamiltonian for a gas with two spin components 



and separable binary interactions d 0, l23il . Consequently, 
we rule out the macroscopic occupation of the Feshbach res- 
onance level as a mechanism behind the superfluid pairing in 
these systems. To this end, we compare predictions obtained 
from both the standard fermionic and boson-fermion Hamil- 
tonians adjusted in such a way that they both recover the low 
energy binary collision physics over a wide range of mag- 
netic field strengths, beyond the regime of universality. These 
adjustments are performed on the basis of a two-channel de- 
scription of the resonance enhanced binary scattering involv- 
ing the entrance channel of the spatially separated atoms and 
a closed channel strongly coupled to it via the Feshbach res- 
onance level. Our approach ll24ll depends on five measurable 
parameters of a Feshbach resonance and is applicable to both 
narrow (closed channel dominated) and broad (entrance chan- 
nel dominated) resonances. 

Throughout this paper we shall discuss a balanced mixture 
of ""'K atoms prepared in the (/ = 9/2, = -9/2) and 
(/ = 9/2, m/ = -7/2) Zeeman states (cf., e.g., Ref. 01). 
We consider magnetic field strengths in the vicinity of the 
202 G (1 G = 10"'* T) zero-energy resonance (singularity of 
the scattering length). This Feshbach resonance is par- 
ticularly well suited to demonstrate the wide range of appli- 
cability of our approach, as its width is in between those of 
the experimentally relevant extremely narrow (closed chan- 
nel dominated) and broad (entrance channel dominated) res- 
onances of ''Li at about 543 lUl and 830 G |7, 8, 1% Q, 
respectively. While the general differences between closed 
and entrance channel dominated Feshbach resonances were 
discussed in Ref. |26|, the very different universal regimes of 
magnetic field strengths of ^Li resonances and their potential 
relevance to the BCS-BEC cross-over were the subject of the 
recent studies of Ref. |27]. 

The paper is organised as follows: In Section|n]we describe 
the two-channel approach to the resonance enhanced low en- 
ergy binary scattering observables and introduce the five rele- 
vant physical parameters for the 202 G zero-energy resonance 
of "^''K. These two-body considerations reveal that the admix- 
ture of the Feshbach resonance level to the experimentally rel- 
evant bound Feshbach molecule never exceeds 8% over the 
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entire range of magnetic fields relevant to the experiment of 
Ref. |9]. Its maximum value is reached at about 7.2 G below 
the position of Bq - 202.1 G 0| of the zero-energy resonance. 
These results are in agreement with exact coupled channels 
calculations |28] from which the two-channel approach was 
originally derived in Ref. i3. Our analysis reveals that a 
proper description of the admixture of the resonance level to 
the multi-channel bound and scattering states requires a two- 
channel approach to explicitly account for the entrance chan- 
nel background scattering potential including its highest ex- 
cited vibrational state. Given the small admixture of the res- 
onance level to the Feshbach molecule, we then determine an 
accurate separable single channel potential, which is suitable 
for the standard fermionic many-body Hamiltonian. This po- 
tential is adjusted in such a way that it accurately describes not 
only the magnetic field dependence of the scattering length but 
also the binding energy of the Feshbach molecule. 

In Section we discuss the implications of the nature of 
the resonance enhanced binary collision physics for the poten- 
tials and parameters of both the standard fermionic and boson- 
fermion Hamiltonians. The procedure of adjusting these ef- 
fective Hamiltonians involves the identification of the two- 
body resonance level with the structureless Bose particles of 
the boson-fermion approach. Its background scattering po- 
tential and inter-channel coupling parameters are then deter- 
mined in such a way that they exactly reproduce the predic- 
tions of the two-channel approach of Ref. 1.241 when applied 
to a pair of atoms. 

Section lTvl brieflv summarises the mean-field approach we 
have used to determine the thermodynamic quantities for both 
the standard fermionic and the boson-fermion Hamiltonians. 
In Section |V] we then apply this mean-field approach to the 
BCS-BEC cross-over in a cold gas of "'"K with a particular 
emphasis on the experimental conditions of Ref. 01 . We show 
that, as expected from our two-body analysis, the standard 
fermionic and boson-fermion Hamiltonians lead to virtually 
the same predictions about the equilibrium physics. In par- 
ticular, our analysis supports the picture of a smooth cross- 
over of the pair size in the entrance channel characterised 
by features similar to those of traditional superconductivity 
flll3,bt.di.y3] and exciton andpolariton ill|23|23.Ell conden- 
sates. We therefore conclude that the occupation of the Fesh- 
bach resonance level is irrelevant to the nature of fermionic 
superfluidity in a cold "^''K gas, in contrast to the findings of 
Refs. I32I l33ll . We predict the position of the zero of the 
chemical potential (cross-over point) at 0.32 G below the zero- 
energy resonance. We have, furthermore, analysed the exper- 
imental pairwise projection technique of Ref. 0] by calcu- 
lating the overlap between the fermionic pairs produced on 
the high field side of the '^''K resonance and the bound Fesh- 
bach molecule on its low field side. Our results indicate that 
fermionic pair condensation phenomena should become invis- 
ible once the gas is prepared at magnetic field strengths further 
than about 0.5 G above the resonance position. Finally, Sec- 
tion lVII summ arises our main conclusions. 



II. MAGNETICALLY TUNABLE INTER-ATOMIC 
INTERACTIONS 

At the low collision energies characteristic of cold dilute 
gases the binary scattering physics can be described by a sin- 
gle parameter of the inter-atomic potential, the s-wave scat- 
tering length a. The experimental technique of Feshbach res- 
onances takes advantage of the Zeeman effect in the atomic 
energy levels to widely tune the scattering length using mag- 
netic fields. In the experiments of Ref. |y] the gas was pre- 
pared as a balanced, incoherent mixture of the lowest ener- 
getic (/ = 9/2, m/ = -9/2) and (/ = 9/2, m/ = -7/2) Zee- 
man levels at magnetic field strengths on the order of 200 G. 
Only pairs of unlike fermions in different Zeeman levels can 
interact via s waves which dominate cold collisions (in the 
absence of resonant phenomena associated with higher partial 
waves). Throughout this paper we shall denote the asymptotic 
5-wave binary scattering channel of a pair of these unlike '^"K 
atoms as the entrance channel, while its associated interac- 
tion potential Vbg will be referred to as the background scat- 
tering potential. As the m f degeneracy of the Zeeman levels 
is removed by a homogeneous magnetic field of strength B, 
the potentials associated with the different scattering channels 
can be shifted with respect to each other. The typically weak 
inter-channel coupling becomes significantly enhanced when 
the magnetic field dependent energy E^^siB) (see Fig.0 of a 
closed channel vibrational state ip^es (the Feshbach resonance 
level) is tuned in the vicinity of the dissociation threshold of 
the entrance channel. Since this threshold coincides with the 
zero of collision energy (zero of energy in Fig.^, its virtual 
degeneracy with £',e,s(B) leads to a resonance enhancement of 
the cold inter-atomic collisions in terms of a zero-energy reso- 
nance in the entrance channel, i.e. a singularity of the scatter- 
ing length at the magnetic field strength Bq described by the 
formula: 

a(B)^a^,\l-j—^J. (1) 

The parameters abg and AB are usually referred to as the 
background scattering length and the resonance width, respec- 
tively. 

The magnetic tuning of the energy of the Feshbach reso- 
nance level distorts not only the scattering continuum but also 
the bound multi-channel molecular energy levels. As illus- 
trated in Fig. [2 the position Bq = 202.1 G of the zero-energy 
resonance, i.e. the magnetic field strength at which the scat- 
tering length has its singularity, coincides with the zero of the 
binding energy Eb{B) of the highest excited multi-channel vi- 
brational state. This Feshbach molecular state (pbiB) persists 
on the low field side of the 202 G zero-energy resonance of 
'^''K where a(B) is positive, while it ceases to exist on the 
high field side of negative scattering lengths. We note that 
the measurable resonance position Z?o is shifted with respect 
to the magnetic field strength B,es (see Fig.^ at which the en- 
ergy £',es(B) crosses the dissociation threshold of the entrance 
channel. In fact, the multi-channel bound state (pbiB) and the 
meta-stable Feshbach resonance level are, in general, rather 
distinct objects, in particular, when the inter-channel coupling 
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FIG. 1: (Color online) Binding energies of the highest excited 
multi-channel vibrational molecular bound states (solid and dashed 
curves), associated with a pair of unlike "'"K atoms in the asymptotic 
(/ = 9/2, m; = -9/2) and (/ = 9/2, = -7/2) Zeeman states, 
versus the magnetic field strength B. The Feshbach molecular state 
(f>b(B) and its energy E^,{B) emerge at the position of the zero-energy 
resonance Bq- The linearly varying energy E^^^{B) of the Feshbach 
resonance level i^res and the energy E-i of the highest excited vibra- 
tional state of the background scattering potential are indicated by the 
dot-dashed and dotted lines, respectively. We note that the measur- 
able resonance position Bo is significantly shifted with respect to the 
magnetic field strength Bres at which the closed-channel resonance 
state energy crosses the dissociation threshold of the entrance chan- 
nel. Here and throughout this paper, we have chosen the zero of the 
energy, at each magnetic field strength B, as the sum of the Zeeman 
energies of a pair of unlike fermions. 



IS Strong. 



A. Two-channel approach 

To accurately describe the resonance enhancement of the 
binary collisions as well as the properties of the Feshbach 
molecule of a pair of unlike '**'K fermions, we employ the 
approach of Ref. 124] , which explicitly includes the entrance 
channel and its strong coupling to a closed channel via the 
Feshbach resonance level. According to the effective treat- 
ment of Ref. 1 24], the general form of a two-channel Hamil- 
tonian matrix 



H2B 



2ch 



W 



(2) 



depends on five physical parameters in addition to the atomic 
mass m: The background scattering potential Vbg of the en- 
trance channel is characterised by the energy E-\ of its highest 
excited vibrational state and the background scattering length 
flbe fsT], while the closed channel potential Va{B) can be ef- 
fectively described by the constant difference dE^^^jdB be- 
tween the magnetic moments of the Feshbach resonance level 
and a pair of asymptotically separated '*'^K fermions in the 



(/ = 9/2, m/ = -9/2) and (/ = 9/2, m/ = -7/2) Zeeman 
states. The inter-channel coupling is accounted for by the off- 
diagonal potential W whose strength and range parameter are 
determined by the resonance width AB and the shift Bo - Bies 
(see Fig.[0. 

At the low relative momenta and associated large de Broglie 
wavelengths characteristic for cold gases, the details of the mi- 
croscopic binary interaction potentials are not resolved. Ref- 
erence 1 24 1, therefore, suggests a convenient effective poten- 
tial matrix that recovers all the relevant low energy physical 
observables of the exact microscopic binary interaction. 



1. Adjustment of the background scattering potential 



Following the treatment of Ref. ||24||, we use the following 
separable form of the background scattering potential: 



Vbg = lrbg)^bgCrbgl ■ 

For the form factor \xbg) we choose the Gaussian ansatv 



(3) 



exp 



(-*) 



(4) 



Here |p) is a plane wave momentum state, normalised as 
exp(/p ■ rlTi)l(2nKfl^. The amplitude ^bg and the range pa- 
rameter (Tbg can be adjusted in such a way that Vbg exactly 
recovers Obg and £_! as follows: The first condition for the ad- 
justments is given by the relationship between the background 
scattering length and the zero-energy T matrix associated with 
the background scattering potential. The exact analytic ex- 
pression for this T matrix given in Ref. f24l] yields: 



flbg 



m ^ 



1 + 



-fbg 



(5) 



The second condition is equivalent to the stationary Schro- 
dinger equation for a bound state of the separable potential 
and reads 12411 : 



^bg [ yfnxe'^' erfcix) - 1 j 



1. 



(6) 



Here erfc(x) - e "~du is the complementary eiTor 

function (sometimes refeiTed to as the Hechenblaikner func- 
tion) with the argument x = V'«|£^-i|o"bg/^- We note that the 
posi tivity of the background scattering length Obg - 174flBohi 
PSl (flBohr = 0.052918 nm is the Bohr radius) guarantees the 
existence of a single bound state of the separable potential Vbg 
for two unlike '^''K fermions. Equations (|5} and (|6j then si- 
multaneously determine the parameters crb„ and ^b^ 



2. Single-resonance approximation in the closed channel 

In most of the experimentally relevant cases the entrance 
channel is strongly coupled only to a single closed channel 
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vibrational state, i.e. the Feshbach resonance level 0res- It is 
therefore usually sufficient to restrict the spatial configuration 
of a closed channel atom pair to (l>r^Ar), where r is the inter- 
atomic separation. To this end Ref. [24] introduces the single- 
resonance ap prox imation (referred to as the pole approxima- 
tion in Ref. 12411 ') which consists in replacing the diagonal 
closed channel part of the general Hamiltonian matrix ^ by 
the one dimensional projection onto the resonance level, i.e. 



(7) 



The resonance energy varies virtually linearly with the mag- 
netic field strength B and can therefore be described by the 
first term of the power series expansion about its zero at Bjes, 
which reads: 



EM 



dB 



'■(B - B,es) . 



(8) 



Here dE^^^ldB is the magnetic moment of the resonance level 
with respect to the sum of magnetic moments of a pair of 
asymptotically separated unlike fermions. It can be obtained 
either from measurements of the binding energies of the Fesh- 
bach molecule away from the zero-energy resonance or from 
exact coupled channels calculations using microscopic inter- 
actions. 



3. Inter-channel coupling 

In accordance with the preceding single-resonance approx- 
imation to the closed channel Hamiltonian the inter-channel 
coupling is determined by the product W|0i-es) |24]. We repre- 
sent this product in terms of an amplitude ^ and a wave func- 
tion i.e. 



W|0,-e,s> = \X)(- 



(9) 



The wave function \x) accounts for the spatial variation of the 
inter-channel coupling in terms of a range parameter cr. We 
choose a Gaussian ansatz which, in the convenient momen- 
tum space representation, is given by: 



exp 



<pUr> = xip) = 



(10) 



Given that the range parameter o-^o of the background scat- 
tering potential has been determined using the procedure from 
above, the parameters ^ and cr can be adjusted in such a way 
that they exactly recover the width AB of the zero-energy res- 
onance as well as its shift Bo - Bres (see Fig.^. Following the 
procedure of Ref. |24], the two conditions for the simultane- 
ous adjustment of ( and cr are given by the exact expressions 

i2 



^B 



47r^2flbg (^^res IdB) 

for the resonance width, and 



flbg 



Bo - Bres = (AB) 



flbg 



sfncr \<tJ 



(11) 



(12) 



for the associated shift. Here we have introduced the mean 
range parameter: 



While AB is usually known from measurements of the 
magnetic field dependence of the scattering length a{B) 
(cf. Ref. 1 6]), the precise value of Bo-Bres can not be easily de- 
duced from experimental data. Reference i24il therefore sug- 
gests an analytic treatment, based on ideas of multi-channel 
quantum defect theory L36t.l , to determine the resonance shift. 
This treatment relates Bq - B^es to the van der Waals dispersion 
coefficient Ce in terms of the mean scattering length 



ispers 



__ 1 r(3/4) l /mC6\'/' 

" V2r(5/4)2l n'- I ■ 



(14) 



Here F is the gamma function. The van der Waals disper- 
sion coefficient characterises the (non-retarded) asymptotic 
behaviour Vbg('") ~ -Ce/r^ of the exact microscopic back- 

r—*oo 

ground scattering potential at large inter-atomic distances r 
and can be deduced from experimental observations 1 38]. Fol- 
lowing the ideas of Ref. |36], the magnitude of Bo - Bres is 
accurately determined by the formula: 



Bo - B,,s = (AB). 



(15) 



The procedure for the determination of the inter-channel cou- 
pling parameters ^ and cr thus consists in calculating the res- 
onance shift using Eq. il5i and then simultaneously solving 
Eqs. O and 



4. Molecular bound states 

The two-channel Hamiltonian of our approach of Ref. f24 
can support two bound states. Their wave functions associated 
with the entrance- and closed-channel components have the 
following general form: 



I Gbg(£b)H'«^ies 



(16) 



Here Gbg(£b) is the energy-dependent Green's function asso- 
ciated with the entrance channel, whose explicit expression 
reads 



Gbg(£b) 



£b- V^ + Vb, 

m 



(17) 



while the normalisation constant Ny, is given by the formula 



Nb= {<Pr,M [Gbg(£b)]' W|0,-es) . (18) 

The energies E\, associated with the two-channel bound states 
in the single resonance approximation are determined by the 
following condition: 



£b = £res(B) + <0i-e,s|WGbg(£b)W|0,-e,s) . 



(19) 
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The energies of the Feshbach molecular state and of the sec- 
ond more tightly bound two-channel state of '^^K are indicated 
in Fig.[nby the solid and dashed curves, respectively. 



5. Parameters for a pair of unlike '^K fermions 

For the experiment of Ref. | 9], which uses a balanced in- 
coherent mixture of '**'K atoms in the (/ - 9/2, = -9/2) 
and (/ - 9/2, nif = -7/2) Zeeman states at magnetic field 
strengths in the vicinity of the 202 G resonance, the phys- 
ical parameters are given by: a^g = 174aBohr |35], /s-i = 
-h X 8.75 MHz 1281, AB = 7.8 G l"^, Bo - Bres = -9.274 
G, and dE.^JdB = 1.679;UBohr (/"Bohr = 9.27400899 x lO^^s 
J/G is the Bohr magneton). We have obtained the reso- 
nance shift Bo - fires from Eq. M5\ using Ce = 3897 a.u. 
f38] (1 a.u.=0.095734 x lO^^^J m'') and the resonance slope 
dEyf-s/dB has been extracted from exact coupled channels cal- 
culations of the binding energy of the Feshbach molecule 
lEsil . From these physical quantities we have determined 
the following parameters of the two-channel Hamiltonian: 
m^bg/(47rft2) - -244.852 flBohr, crbg = 57.387 aBohr, cr - 
33.123 flBohr, and m^^ HAnh^a-) = /i x 80.272 MHz. Figure 
[^shows the binding energies of the multi-channel vibrational 
states predicted by the two-channel Hamiltonian of Eq. us- 
ing the parameters listed above. We note that the two-channel 
approach also exactly recovers the magnetic field dependence 
of the scattering length a{B) described by Eq. Q. 



B. Single-channel approach 

As the magnetic field strength approaches the position Bq of 
the zero-energy resonance from the positive scattering length 
side, the admixture of the resonance level 0ies to the physi cally 
relevant Feshbach molecular state 4>h{B) vanishes I 22l l2al . 
In this universal regime of magnetic field strengths it is al- 
ways possible to introduce an effective single-channel Hamil- 
tonian which accurately recovers the resonance 
enhancement of the scattering as well as the universal binding 
energy 



(20) 



of the Feshbach molecule. The studies presented in Refs. f2l', 
I26I1 have shown, however, that the range of validity of such 
a single-channel approach can be significantly extended by 
properly adjusting the binary potential of the effective Hamil- 
tonian, provided that the resonance is broad and entrance 
channel dominated iijll . We shall show below that the 202 G 
resonance of '^'^K fulfils all requirements for an accurate treat- 
ment of the binary low energy bound and continuum spectrum 
in terms of just a single asymptotic scattering channel. 

The single-channel approach of Refs. |2ll 13 l2^ eff'ec- 
tively describes the resonance enhancement of the low energy 
scattering as a perturbation of the background scattering po- 
tential. Following this approach, we replace the parameter ^bg 



in Eq. (|3} by a magnetic field dependent amplitude ^eff-ich(fi) 
and insert the associated eff'ective single-channel potential 

Veff-Ich(B) = lrbg>^eff-Ich(B)(Afbgl (21) 

into the general form of a single-channel Hamiltonian 

i/2B-lch = V2 + yeff-lch(B) . (22) 

m 

In this treatment the range parameter cr[,g retains the same 
value as in the two-channel approach, while the amplitude 
^eff-ich(B) is adjusted in such a way that the eff'ective poten- 
tial recovers, at each value of the magnetic field, the scatter- 
ing length a{E) of Eq. Q. This adjustment can be performed 
in complete analogy to the determination of ^bg in the two- 
channel approach of Ref. fl^. The counterpart of Eq. (|5} for 
the determination of ^eff-ich(fi) leads to the relationship: 



^eff-lch(fi) - 



4nfp-a(B)/m 

1 - fl(B)/( V^CTbg) 



(23) 



We note that despite the divergence of the scattering length 
a(B) at resonance the amplitude of the effective single-channel 
potential ^eif-ich(fi) remains a smooth function of the mag- 
netic field strength B. The single-channel approach not only 
recovers the exact magnetic field dependence of the scatter- 
ing length, by identifying a{B) in Eq. ( I23> with Eq. Q, but it 
also leads to the exact asymptotic limit E^i of the binding en- 
ergy £'b(B) away from the resonance (see Fig.[0, by using the 
same form factor \xbg} as the separable background scattering 
potential of Eq. Q. 

Figure 12 compares the binding energy £'b(B) of the Fesh- 
bach molecule, as obtained from the two-channel approach of 
Ref. ll24ll . with the predictions of the effective single-channel 
approach, as well as analytic estimates. These estimates are 
based on Eq. (|5DJ and on the formula 



m(a — a)^ ' 



(24) 



which follows from the treatment of Ref. ^yh and accounts 
for the non-retarded van der Waals tail ybg(r) ~ -Ce/r*' of 

the background scattering potential, in terms of the mean scat- 
tering length a of Eq. MA\ . in addition to the magnetic field 
dependence of the scattering length a(B). As the magnetic 
field strength B approaches the position Bq of the zero-energy 
resonance from below, all predictions coincide. In this uni- 
versal regime of magnetic field strengths the scattering length 
very much exceeds all the other length scales set by the binary 
interactions and completely determines the magnitude of the 
binding energy. Further away from the resonance the qual- 
ity of the universal estimate of Eq. ( I20> deteriorates, while 
Eq. MAI provides a good approximation of Et,{B) over the en- 
tire range of magnetic field strengths shown in Fig.|2 The full 
predictions of the two- and single-channel approaches are vir- 
tually indistinguishable. This agreement suggests that the ad- 
mixture of the resonance level 0res to the Feshbach molecule 
remains small also outside the universal regime of magnetic 
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FIG. 2: (Color online) Binding energy E^,{B) of the Feshbach molec- 
ular bound state 0b(B) associated with a pair of unlike ^''K atoms 
in the asymptotic (/ = 9/2, m/- = -9/2) and (/ = 9/2, ot/ = -7/2) 
Zeeman states versus the magnetic field strength B. The figure shows 
a comparison between the results of the two-channel approach (solid 
curve), the effective single-channel approach (dashed curve), and an- 
alytic estimates of the binding energy (dotted and dot-dashed curves). 
While the dotted curve refers to the universal estimate of Eq. <20t . 
just depending on the scattering length a(B), the dot-dashed curve 
also accounts for the non-retarded van der Waals interaction -Ce/f* 
of the background scattering potential, at asymptotically large inter- 
atomic distances r, in terms of the mean scattering length a I2lll37ll 
via Eq. <24l . 



field strengths, i.e. the 202 G resonance of '^''K is entrance 
channel dominated. In fact, we shall show in Section fVl that 
the closed-channel admixture to the bound state wave func- 
tion 4>h{B) never exceeds about 8% and reaches its maximum 
at about 7.2 G below the position Bo - 202.1 G of the zero- 
energy resonance. 



III. MANY-BODY HAMILTONIANS 

A. The microscopic Hamiltonian 

Having established the microscopic parameters of the inter- 
atomic interactions, which properly describe the experimen- 
tally relevant part of the two-body energy spectrum, we shall 
now discuss the many-body Hamiltonians that we use in our 
studies of the BCS-BEC cross-over in a cold dilute gas of 
"*"K. The most general many-body Hamiltonian, including all 
atomic internal degrees of freedom, has the following form: 



ek,v|flk^, flk,, 



vi,k 



^Vi,V2,V3,V4(k, k ) 



2 

X d 



a+k.vi'^a-k.v,"?- 



k',V3"f +k',V4 • 

(25) 



Here the indices v\, V2, V3, V4 label the atomic Zeeman lev- 
els. The wave vector k in the kinetic energy contribution to 
Eq. (I25> is associated with a single-particle momentum, while 
the pairs of wave vectors q and k in the interaction term corre- 
spond to two-body centre-of-mass and relative momenta. All 
these indices refer to our choice of plane wave basis states 
exp(/k ■ r) / in a box of volume 'V with periodic boundary 
conditions. The associated annihilation and creation operators 
obey the usual fermionic anti-commutation relations: 



-<Jv,,v,<Jk,k' 



=0. 



(26) 
(27) 
m Eq. M5\ 



«k,v, + «k',V2"k,, 
ak,v,ak',V2 + ak',V2flk,)_ 

We note that the microscopic potentials Vvi.vt.vj.vj 
all decay to zero in the limit of large inter-atomic separation. 
The atomic Zeeman energies, which determine the dissocia- 
tion thresholds of the different two-body asymptotic scattering 
channels are included in the single-particle energies ek.v, ■ 

B. Ambiguities between diflferent many-body approaches 

While a mean-field treatment of the BCS-BEC cross-over 
using realistic microscopic interactions for all binary scatter- 
ing channels seems impractical, there are several ways to for- 
mulate approximate theories, derived from the Hamiltonian of 
Eq. M5\ or the effective boson-fermion model flTlfT^l . which 
are all compatible with the two-channel approach to the binary 
collision physics of Section |ll] One possible condition for 
such a compatibility may be, for instance, that the binding en- 
ergy of the Feshbach molecule, predicted by the two-channel 
approach, exactly coincides with twice the chemical potential 
obtained from the mean-field theory in the asymptotic limit of 
zero density. 

During the course of our studies we have derived such a 
mean-field approach 139], based on the microscopic Hamilto- 
nian of Eq. \25\ and the general kinetic theory of Refs. |43 
I41I1 . by applying the two-channel approach of Section [IIJ in- 
cluding the single-resonance approximation, to the two-body 
propagators. This treatment is most sensible as it relies upon 
controlled approximations rather than uncontrollable model 
assumptions about the form of an effective Hamiltonian. The 
dynamic equations for the one-body density matrix, however, 
explicitly depend on the single-atom Zeeman states that con- 
stitute the closed channel associated with the Feshbach reso- 
nance level. The atomic energy levels need to be explicitly 
specified in such a description once one and the same single- 
particle Zeeman state is shared between the relevant two-body 
entrance and closed channels. This phenomenon has been 
discussed previously in Ref. |42] for the case of '^"K and re- 
flects physical intuition, since measurements of single-particle 
quantities, like the position of an atom or its internal state, do 
not necessarily provide information about whether this atom 
was coiTelated with another partner in the spatial configura- 
tion of the Feshbach resonance level. 

The assumption of such a generally unphysical possibil- 
ity of distinguishing spatially separated from correlated atom 
pairs in a gas on the basis of one-body observables under- 
lies the boson-fermion model 1 17nl8>l . This model separates 
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out the Feshbach resonance state on the level of the Hamil- 
tonian. We shall show, however, that the standard thermody- 
namic physical quantities we discuss in Section |y] are recov- 
ered, at least on the level of the mean-field approximation, 
by many-body approaches based on both the boson-fermion 
and the standard fermionic Hamiltonians, provided that they 
are suitably adjusted to be compatible with the precise two- 
channel approach to the two-body physics of SectionUl] This 
observation suggests that the thermal equilibrium properties 
of Section |V] are insensitive to the underlying model Hamil- 
tonian, regardless of the level of accuracy to which they in- 
clude the separation into different asymptotic binary scatter- 
ing channels. 

C. Adjustment of the boson-fermion model 

The boson-fermion model leads to an effective two-channel 
description of the two-body physics and includes only the in- 
teractions between unlike fermions in the different Zeeman 
states associated with the entrance channel of Section|ll] The 
boson-fermion Hamiltonian has the following form: 

i,k q 

+ ^bg('^''*^')«2+k,T^l-k,i"l-k'.i"!+k'.T 
k,k',q 2 • 2 ■ 

+ ^ §(k) (^;,fl|_k,i«f +k,T + /'•c-) ■ (28) 

k,q 

Here and b\^ annihilate or create pairs of atoms with the 
centre-of-mass momentum q in the Feshbach resonance level, 
respectively, and fulfil the bosonic commutation relations: 

- ^q'^q = ^qq' ' (^9) 

- K'K = . (30) 

The kinetic energy associated with the centre of mass of an 
atom pair in the meta-stable Feshbach resonance level is de- 
noted by Eq. In analogy to related applications in the theory 
of superconductivity, we choose, in the following, the notation 
s -1, i to label the single-particle Zeeman levels that consti- 
tute the 5-wave entrance channel. We assume all fermionic 
operators Ok.j and ^ to commute with all bosonic operators 

bq and £q. In analogy to our considerations of the two-body 
problem in Section|ll]we have included the magnetic field de- 
pendent single-particle Zeeman energies in the binary interac- 
tions. The single-particle energies ek are therefore indepen- 
dent of the index s associated with atomic Zeeman levels. 

In order to determine the parameters of the boson-fermion 
Hamiltonian, we compare the Schrodinger equations obtained 
from Eq. ( I28> and the binary two-channel Hamiltonian of 
Eq. (0 when applied to a general two-body state with com- 
ponents in both the entrance channel and the Feshbach reso- 
nance level. As both approaches are required to describe the 
same physics, this comparison shows that £'ies(B) is the res- 
onance energy given by Eq. (|8}, while Vbg is the background 



scattering potential of Eq. (|3}- Its matrix element in terms of 
the energy states of the periodic box is thus given by 

Vbg(k,k') = f dr f flfr'e-'''--ybg(r,r')e'''' '"', (31) 

y j-v J-v 

where the integration extends over the volume of the box. 
The parameter 

g(k) = -L r flfre-'"--W(r)0res(r), (32) 

which characterises the inter-channel coupling, can be ob- 
tained from Eq. We note that with these adjustments of 
potentials and parameters the boson-fermion Hamiltonian of 
Eq. ( I28> and the two-channel approach in the single resonance 
approximation of SectionlUprovide exactly identical descrip- 
tions of the two-body physics. 

D. Adjustment of the standard fermionic Hamiltonian 

A many-body approach compatible with the single-channel 
description of resonance enhancement in the binary low en- 
ergy collision physics of Section HI] can be obtained from a 
many-body Hamiltonian usually applied in the theory of su- 
perconductivity ,4^ .23J . This standard fermionic Hamilto- 
nian is given by the following expression: 

^MB-F = ^fkak,A.* + ^eif-lch(k,k',B) 
sM k,k',q 

^ ^l+k.T^l-k.i^i-k'.i^f +1''.T ■ ^^^^ 

Here Veff-ich(fi) is the interaction potential for a pair of atoms 
in the Zeeman states associated with the entrance channel, 
which we have labelled by the indices s =1,1- We choose 
this potential to be identical to the binary interaction in the 
single-channel Hamiltonian of Eq. ( I22t . This choice assures 
that the two-body physics described by the standard fermionic 
Hamiltonian exactly agrees with the single-channel approach 
of Sectionin 

As both the two- and single-channel approaches of Section 
m] yield virtually the same low energy two-body physics, we 
may also expect a similar agreement between the predictions 
obtained from the Hamiltonians of Eqs. J28> and J33> with 
regard to the thermal equilibrium physical quantities of Sec- 
tion|V] The formal establishment of this equivalence of these 
many-body Hamiltonians in applications to standard mean- 
field theory will be the subject of the following sections. Since 
we have included smooth binary potentials with a proper spa- 
tial extent in the Hamiltonians, our mean-field calculations do 
not involve ultra-violet divergences and, therefore, do not re- 
quire any renormahsation procedures. 

IV. MEAN-FIELD APPROACH TO THE THERMAL 
EQUILIBRIUM 

At temperatures close to absolute zero or well below the 
critical temperature T^, the BCS-BEC cross-over has been 
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successfully analysed by mean-field theory in terms of the 
BCS wave function, which smoothly interpolates between the 
BCS and BEC regimes flUJ^. The derivation of the ther- 



mal equilibrium properties [32 3 3', "43*, "44", 45 , 4?, "i?, "i?, "i?, 
EO. IL 52. 53. 54, 55, 56, 52, ,58, 59, 60, 6 L 62, 63. M. 65J 
follows standard techniques (see, e.g., Ref. 12311 '). We thus in- 
clude just a brief derivation of their extension to the boson- 
fermion model, applying path integral methods |66]. The 
associated mean-field approach for the standard fermionic 
Hamiltonian model can be obtained from our analysis simply 
by replacing g{k) — > and Vbg — > Veff-ich in all formulae. 



A. Path integral approach 

All thermodynamic equilibrium physical quantities of a gas 
can be derived from the quantum partition function, which is 
given by the general formula: 



Z = ti- [^ 



,-/3(Hmb-bf-/jW)1 



(34) 



Here the factor /3 - llik^T) accounts for the temperature T 
scaled to energy units by the Boltzmann constant ^b. ^ is the 
number operator and ji the chemical potential, while "tr" indi- 
cates the trace of an operator. The quantum partition function 
Z can be conveniently represented in terms of a coherent state 
path integral. To this end, we introduce the action 



S = 



dT 



(35) 

where is the total number of atoms and //mb-bf is the many- 
body boson-fermion Hamiltonian of Eq. j28> . In this Hamilto- 
nian flk.i and bq are replaced by r-dependent Grassmann num- 
bers flk.s and c-numbers bq, respectively. We have, further- 
more, introduced the quantities flk,; 
spinors 



bq as well as the Nambu 



(36) 



in Eq. J35L where the bar indicates that the quantities flk.s 
and flk.s. etc., are independently varied. In our considerations 
flk.s and flk.s as well as bq and bq are all related by complex 
conjugation. The coherent state path integral then determines 
quantum partition function Z in terms of the action to be: 



z 



, i/r, b, b)e' 



(37) 



B. 



Relationship between the path integral approach and the 
zero-temperature mean-field approximation 



In the following, we shall briefly summarise the assump- 
tions underlying our evaluation of Eq. i'ili in the mean- 



field approximation as well as their relationship to derivations 
based on the condensate wave function of Refs. LI„^|5|]. In 
order to treat the interaction we isolate the components as- 
sociated with the pairing. Such a procedure is equivalent to 
considering only those contributions of the total Hamiltonian 
that yield non-zero matrix elements in the condensate state, 
i.e., (OI^mb-bfI'I') ?^ 0. Following the ideas of standard BCS 
theory |23], the associated zero-temperature ground state con- 
densate wave function is given by the formula: 



exp (-1 AP/2 + Abbl + 2k WkSk T^-k i) 
I*) = ■ '■ ^|vac) 



nk'(i + iwk'p) 



2^l/2 



)|vac). (38) 



Here |vac) is the vacuum state. We have, furthermore, intro- 
duced the quantity = Vk/Mt in Eq. j38> which depends on 
the quasi-particle amplitudes Mk and Vk- The quasi-particle 
amplitudes fulfil the requirement |Mkl^ + |vkP - 1, where 

2 Ivkl' = <1>l «^ak,T + al^a^,i) |0) (39) 

determines the average occupation number of the momen- 
tum mode associated with the wave vector k. The prod- 
uct % = Mki'k can be interpreted as a pair wave function, 
whose modulus squared, summed over all wave vectors k 
(i.e. 2k Iwkt^'kP), provides a measure for the number of con- 
densed fermionic pairs. Within the framework of the boson- 
fermion model, the coefficient A\, in Eq. ( I38> is related to the 
average number of pairs of atoms in the resonance state con- 
figuration through the formula: 



(40) 



We shall omit, in our evaluation of Eq. J37> . the decou- 
pling in the diagonal, Hartree-Fock, channel because, firstly, 
the Fock pairing, which results from exchange, requires indis- 
tinguishable atoms in the same Zeeman states. The associated 
binary interactions, however, which all involve non-isotropic 
partial waves, are omitted from the outset in the i-wave boson- 
fermion Hamiltonian. Secondly, the inclusion of the Hartree 
term results in only a negligible energy shift. 



C. Hubbard-Stratonovich transformation 



To further evaluate the finite temperature partition func- 
tion of Eq. ( I37> . we proceed by applying the Hubbard- 
Stratonovich decoupling to the interaction in the off-diagonal 
channel, which yields: 
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(41) 



r 



The transformed partition function depends on the or- 
der parameter Ak, sometimes referred to as the Hubbard- 
Stratonovich field, whose explicit expression within the mean- 
field approximation in the limit of zero temperature reads 
= -2k' ^bg(k,k')<<l>|fl-k',iflk',Tl4>)- In our more gen- 
eral finite temperature treatment the order parameter is a t- 
dependent function which we shall determine using a saddle- 
point analysis. We have, furthermore, introduced in Eq. (I41> 
the quantity 



rP 

5o= flfT2A*ybg'(k,k')Ak, 

^0 k,k' 

Jo 



(42) 



where V^^ indicates the inverse potential energy operator as- 
sociated with the background scattering in the relative mo- 
tion of a pair of unlike fermions. We note that, in the mean- 
field approximation, the partition function of Eq. \A\\ depends 
only on the c-numbers feq=o and l7q=o associated with the zero- 
momentum mode of the centre-of-mass motion of a pair of 
fermions in the resonance level. 

Since the exponent in Eq. (I41> is a quadratic form in the 
fermionic fields, the associated Gaussian integral can be per- 
formed straightforwardly. This integration leads to the follow- 
ing formula: 



z 



This formula depends on the thermal Green's function G, 
whose inverse operator is given by: 



G ' =5^ - (fk - /^) 0-3 - [g{y.)bo + Ak] 0-+ 
- [§(k)/7o + Ak] cr_ . 



(44) 



Here cr\ , 0-2 and 0-3 are the usual two-dimensional Pauli spin 
matrices, i.e. their commutation relations are given by cr\cr2 - 
(Ttcri - 2i(j-i and its cyclic permutations, while cr+ = (cri + 
io-2)/2 and cr_ = (en - io-2)/2 are the associated raising and 
lowering matrices. 



D. Saddle-point analysis 

We shall apply, in the following, a saddle-point analysis in 
which the order parameters Ak(r) and boir) are assumed to 
be independent of t. This approach minimises the free en- 
ergy as the derivatives of Ak and bo, with respect to t, do not 
contribute. The Green's function G is then diagonal in the 
Matsubara frequencies and in the wave vectors k. We shall. 



therefore, employ a discrete Fourier analysis to change the 
representation from the variable r to the Matsubara frequen- 
cies (jjj = {2j + \)nlfi, where j is an integer Performing the 
Fourier expansion of G after inversion, the mean-field ther- 
mal Green's function is given by the formula: 

Gkojj = A - iojj - (fk - y^)o'3 - [g(k)^o + Ak] 0-+ 

-[g(k)bQ + Al]cr_}. (45) 



Here Zsk - ^J(ek - jJ^)^ + \g(}Obo + AkP are the single-particle 
excitation energies. Varying the action with respect to A^ 
and bo leads to the saddle-point (mean-field) equations. We 
then perform the standard Matsubara frequencies summation 
and consider the uniform gas limit, obtained by the follow- 
ing substitution Ek ^ / '^(^ - f^J^^- I" "^^e follow- 
ing, we shall use the explicit form of the potentials of Sec- 
tion HH i.e. the background scattering potential Vbo(p,p') = 
^hgXhgip)Xbgip') and the inter-channel coupling g{p) = ^^(p), 
in addition to the abbreviation 



dpXbgip)Ap = ^bgA . 



(46) 



The saddle-point equations are then given by the following 
pair of formulae: 



= jD(A*,A,^o,^o)e"^°exp(tr[lnG"']) . (43) 



[£,es(B) - 2iu] bo = jdp 



2V(ep-//)2 + m;,)P \ 2 



(47) 



Here the quantity 

np) ^ ^X(p)bo + ^bg)(hgip)A 



(48) 



(49) 



may be interpreted as a self energy. The chemical potential ju 
is determined by the density equation 



1 



(2nny 



2bobo + j 



dp 



1 - 



^|W^WTWpW) 



(50) 

where n is the number of atoms per unit volume. The mean- 
field partition function reads: 



ZMF = e-^°exp(tr[lnG-^]) . 



(51) 



Here So and G given by Eqs. (I42> and i44\ . respectively, 
are now evaluated at the saddle point, i.e. bo and A are the 
solutions of Eqs. (|47j, (|48j and 
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E. Solutions of the saddle-point equations 

In SectionlVlwe provide the numerical solutions of the three 
coupled equations ( I47> . ( I48> and (15 0> . To obtain a qualitative 
overview of their behaviour, we note that in the zero-density 
limit n — > the chemical potential yU recovers half the molec- 
ular binding energy and the equations i47i . ( I48t and ( 15 Ot re- 
duce to the two-body Schrodinger equation associated with 
the two-channel Hamiltonian (|5}- There are, however, two 
different solutions to this equation. One of them corresponds 
to the weakly bound Feshbach molecular state just below the 
entrance channel dissociation threshold (solid curve in Fig- 
ure^. This state can be efficiently populated by an adiabatic 
sweep of the magnetic field strength. The second solution cor- 
responds to the next, more tightly bound, molecular state of 
energy given by the dashed curve in Figure^ Similarly, at the 
many-body level (i.e. n + Q) there are three different extrema 
of the effective action in Eq. ( I43> . which are all solutions of 
Eqs. ( I47L ( I48> and ( 15 0> associated with different chemical po- 
tentials. The trivial uncondensed solution (A - Q,b(, - 0) is 
unstable below T^, but there are two stable solutions with a 
finite condensate component. One of them corresponds to the 
condensation of weakly bound Feshbach molecules with an 
energy just below the entrance channel dissociation threshold 
(solid curve in Figure^. The other solution corresponds to 
the condensation of molecules in the more tightly bound state 
(dashed curve in Figure [O- The latter solution is given by a 
lower minimum of the effective action in Eq. ( I43> . but it is not 
populated on the time scale of the experiments of Ref. ||9j|. In- 
stead, the local minimum (quasi-equilibrium) solution is the 
physical one, relevant to the experiments of Ref. [gj. We shall 
analyse this solution in SectionfVl 

The saddle-point equations for the mean-field approach as- 
sociated with the standard fermionic Hamiltonian of Eq. (I33> 
can be easily obtained from Eqs. (I47> . ( I48> and (I50> in the 



limit g(p) — > 0, in addition to the replacement Vbg — > Veft-ich- o 



The single-channel effective potential Veff-ich supports just 
a single bound state, the Feshbach molecule, which is effi- 
ciently populated in the experiments of Ref. |9]. The associ- 
ated mean-field equations, therefore, have only two solutions: 
the uncondensed one, unstable below T^, and the condensate 
of weakly bound molecules, which we shall analyse in Section 
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V. BCS-BEC CROSS-OVER IN A DILUTE GAS 



Using the microscopic parameters determined in Sectionllll 
we have numerically solved the coupled saddle-point equa- 
tions ( I47> . ( I48l l and ( 15 0> for a wide range of magnetic fields 
and atomic densities. In particular, we have considered mag- 
netic field strengths from 1 G above to 10 G below the po- 
sition of the zero-energy resonance (Bq - 202.1 G |9]) in a 
mixture of fermionic "*"K atoms in (/ = 9/2, mj = -9/2) 
and (/ - 9/2, my = -7/2) Zeeman states. This range covers 
all magnetic field strengths relevant to Ref. |9] and includes 
both the universal regime around the zero-energy resonance 
(extending to about 1 G below Bq) and the region in which 



the energy of the Feshbach molecule is strongly influenced 
by the anti-crossing between the energies of the two highest 
excited diatomic multi-channel vibrational bound states (solid 
and dashed curves in Fig.^. We have also considered a wide 
range of experimentally accessible densities from 10'^ cm"^ 
to 5 X 10'"* cm"^. In order to establish a characteristic value of 
the density of atoms in a homogeneous two-component gas, 
we require it to reproduce the Fermi energy of the trapped gas 
for the conditions of Ref. |9]. This characteristic density is 
equal to 1.5 x lO'-' cm"-' and is given by the two-component 
peak density of atoms in the trap. 



A. Standard fermionic versus boson-fermion approach 

We have performed the numerical analysis for both ap- 
proaches, using the standard fermionic Hamiltonian and the 
boson-fermion Hamiltonian, and found that in the entire range 
of experimentally relevant magnetic field strengths and densi- 
ties there is hardly any difference visible between them for 
the standard many-body observables considered in this pa- 
per Consequently, in the following discussion, we just give 
the results obtained from the standard fermionic Hamiltonian 
and explicitly point out small differences with respect to the 
boson-fermion approach when they arise. 




FIG. 3: (Color online) The fraction of the population associated with 
the Feshbach resonance level of the closed channel as a function of 
the magnetic field strength B and the density n. The solid curve cor- 
responds to the density characteristic for the experiment of Ref. |^ 
while the result for the lowest density studied (the dot-dashed curve) 
is practically indistinguishable from the closed-channel admixture 
to the coupled-channels bound state wave function of the two-body 
problem. The inset enlarges the region in the vicinity of the zero- 
energy resonance. We note that, regardless of the density, the maxi- 
mum population associated with the bosonic field operator bo is only 
8% for B - B(, ~ -7.2 G. Given the virtual density independence of 
this maximum population, this result also coincides with the maxi- 
mum closed-channel admixture to the highest-excited vibrational di- 
atomic bound state of the two-channel two-body Hamiltonian J2j. 

At the two-body level, the similarity of the results obtained 
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from both many-body approaches is principally due to the 
fact that the admixture of the closed channel configuration 
to the highest excited multi-channel vibrational bound state 
is rather small for the entire range of experimentally relevant 
magnetic field strengths. In the many-body analysis of the 
boson-fermion approach this phenomenon is reflected in the 
small population associated with the bosonic field operator 
bo, as shown in Figure |3 In the vicinity of the position of 
the zero-energy resonance (i.e. for magnetic field strength dif- 
ferences B - Bo from -1 G to 1 G) this population exhibits 
some density dependence, but its magnitude is always less 
than 4%. In particular, for the density relevant to the experi- 
ment of Ref. |i9|l the Feshbach resonance population above Bq 
reaches at most 0.6%, which is about one order of magnitude 
lower than the predictions of Refs. |32, 33]. Further below 
Bq, it becomes essentially density independent and reaches a 
maximum of about 8% at 7.2 G below the zero-energy res- 
onance. This result agrees with the predictions of two-body 
coupled-channels calculations |28|. For even smaller values 
of the magnetic field strength the population associated with 
bo decreases. At the two-body level, this phenomenon stems 
from the influence of the anti-crossing between the two high- 
est excited vibrational bound states (solid and dashed curves 
in Figure 0. The domination of the entrance channel com- 
ponents in the states of weakly-bound molecules reflects their 
long-range nature 1 22|] . 



B. Pair wave function and single-particle excitation spectra 

Several theoretical studies d H, |^ have predicted that 
the cross-over between the BCS phase of correlated fermion 
pairs and the BEC phase of tightly bound diatomic molecules 
may be realised either by increasing the strength of the 
inter-particle interactions (the route employed in dilute Fermi 
gases) or by decreasing the particle density (as in the case 
of exciton and polariton BECs). In the following, we shall 
illustrate, as a function of density and interaction strength, 
the different pairing phenomena in terms of thermal equilib- 
rium physical quantities that can be derived from the mean 
field approach of Section Hvl We shall show, in particular, 
that the BCS-BEC cross-over point (zero of the chemical po- 
tential) in the experiments of Ref. |9] is shifted with respect 
to the position of the zero-energy resonance (singularity of 
the two-body scattering length) towards lower magnetic field 
strengths. This observation corroborates the conventional 
character of the BCS-BEC cross-over in "^''K as op pose d to 
the ideas of resonance superfluidity outlined in Refs. 133.13311 . 

The equilibrium quantities that we shall consider involve 
the pair wave function 



2(k) 



2V(ek-yu)2 + mk)p 



the fermionic distribution function 



1 - 



V(ek-/^)2 + mk)p 



(52) 



(53) 



and the single-particle excitation spectrum 

£k= 7(6k-/i)2 + mk)|2. 



(54) 
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FIG. 4: (Color online) The single-particle excitation spectrum (the 
upper panel), given by Eq. <54> . and the fermionic momentum dis- 
tribution function (the lower panel), determined by Eq. <53t . for the 
density n = 1.5 x 10''' cm"^ and a range of magnetic field strengths 
B, relevant to the experiment of Ref. [9]. We note that the minimum 
of the excitation spectrum evolves from its position ?Llk = on the 
BCS side to = at B - Bo = -0.32 G. 

The single-particle excitation spectra for densities of the 
experiment of Ref. |9] and different values of the magnetic 
field strength (the upper panel of Figure]^ show qualitatively 
different forms in the BEC and BCS limits of the cross-over 
problem. On the BCS side the excitation spectrum has a min- 
imum at non-zero k - (the two lower curves in the upper 
panel of Fig.|3}, where is the Fermi wave number. On the 
BEC side (the two upper curves in the upper panel of Fig.O 
iik recovers the binding energy of a molecule at its minimum 
at A; = 0. For finite wave numbers k the spectrum has a 
quadratic form associated with the dispersion relation of the 
kinetic energy of the molecules. We note that at a finite den- 
sity the position of the zero-energy resonance Bo is not a char- 
acteristic field strength with respect to the cross-over problem. 
In fact, at B = Bo the spectrum shows BCS-like features due 
to its minimum at non-zero k. The value of the magnetic field 
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strength where the minimum of the spectrum changes from 
non-zero to zero k depends on the density. In particular, for 
the density of the experiment of Ref. |9] this characteristic 
magnetic field strength is located at approximately 0.32 G be- 
low Bq. The fermionic momentum distribution function |vkP, 
as given by Eq. i53\ . is shown in the lower panel of Fig.|4] It 
smoothly evolves, as a function of the magnetic field strength, 
from the step-like form of a weakly-interacting Fermi gas (the 
BCS phase) towards a flat distribution characteristic for the 
macroscopic occupation of the lowest energy mode in a con- 
densed Bose gas (the EEC phase). The qualitative change of 
pairing phenomena between the BCS and EEC limits is most 
intuitively reflected in the spatial form of the pair wave func- 
tion. Its oscillatory behaviour on the BCS side is characteris- 
tic for unbound correlated pairs, while its exponential decay 
in the BEC limit recovers the long range asymptotic form of 
the molecular bound state wave function (see Fig.|5|l. 
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FIG. 5: (Color online) The spatial form of the pair wave function 
for the density n = 1.5 x 10'^^ cm"^ and a range of magnetic field 
strengths B. We note that for B - So smaller than about -0.3 G the 
wave function loses its oscillatory form characteristic for the Cooper 
pairs and decays as a function of the inter-atomic distance. Such 
a virtually exponential decay is characteristic for diatomic bound 
states. 

In the case of contact interactions, which are determined 
by a single parameter involving the scattering length, Leggett 
has shown fy] that the qualitative changes in the form of the 
pair wave function and the excitation spectrum occur when the 
chemical potential jj crosses zero. The chemical potential is 
positive on the BCS side and negative on the BEC side of the 
cross-over. In the case of spatially extended or non-local po- 
tentials (such as separable interactions) the qualitative change 
in the form of the pair wave function and the excitation spectra 
does not, in general, coincide with the vanishing of the chem- 
ical potential. Our results indicate, however, that in a dilute 
gas of '**^K atoms at typical experimentally accessible densi- 
ties this boundary practically coincides with Leggett's univer- 
sal boundary at = 0. It is therefore instructive to study the 
chemical potential as a function of the magnetic field strength 
and as a function of the density (see Figure|6}. Far below Bq 
the chemical potential is negative, independent of the density 



and its limiting value recovers one half of the energy of the 
highest excited vibrational molecular bound state. The mag- 
netic field strength where the chemical potential crosses zero 
(the cross-over point) depends on the density. Given the den- 
sity n = 1.5 X 10'^ cm"-', characteristic for the experiments of 
Ref. the cross-over point is located at 0.32 G below the 
position of the zero-energy resonance. 
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FIG. 6: (Color online) The chemical potential /j as a function of the 
magnetic field strength B for different densities n. Away from So, 
in the BEC limit, 2/j approaches the energy of the highest excited 
vibrational molecular bound state £b, as depicted in Figure |2| We 
note that for the density n = 1.5 x 10'^ cm"^ of the experiment |3l 
the chemical potential crosses zero at S - So = -0.32 G. The inset 
enlarges the region in the vicinity of the zero-energy resonance. 

To interpret the results of the pairwise projection technique 
employed in Ref. f^, we have also analysed the magnetic field 
and density dependences of the coherence length ^. This typ- 
ical length scale characterises the size of the fermionic pairs 
and is given by: 

f = —r T- . (55) 

/ dr [/f(r)]2 

Figure0shows the coherence length of the fermion pairs and 
the magnetic field strength for which the size of the pairs 
equals the mean inter-atomic spacing. For the characteristic 
density of the experiment of Ref. |9] the size of the fermion 
pairs equals the mean inter-atomic spacing at a magnetic field 
strength of about 0.5 G above the position of the zero-energy 
resonance. This coincides with the boundary value of the 
magnetic field strength at which condensation phenomena 
were observed via the pairwise projection technique employed 
in the experiments of Ref. iQl • 



C. Pairwise projection of fermionic condensates onto 
molecules 

In the experiments of Ref. ^ a new method of probing 
the paired fermionic component of the gas was used. This 
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FIG. 7: (Color online) The coherence length ^, characterising the size 
of the fermion pairs, as a function of the magnetic field strength B for 
different densities n. The circles indicate the values of the magnetic 
field strength at which the size of the pairs equals the mean inter- 
atomic spacing. We note that for the density n = 1.5 x 10" cm""* of 
the experiment of Ref. |9] (solid curve) the size of the pairs equals 
the mean inter-atomic spacing at B ^ Bq + 0.5 G. 



technique involved a fast ramp of the magnetic field strength 
across the 202 G resonance of '**'K. Ideally, this rapid pro- 
cess projects the centre-of-mass momentum distribution of 
the correlated fermionic pairs on the high field side of Bq 
onto the centre-of-mass momentum distribution of weakly 
bound molecules on the low-field side of the zero-energy reso- 
nance L671 . Such a trapped fermionic condensate of correlated 
pairs is characterised by a narrow centre-of-mass momentum 
spread, which was indeed observed as a distinct peak in the 
measured molecular momentum distribution |9j. In the prac- 
tical implementation of this procedure reported in Ref. 01, 
the magnetic field strength was swept to about 10 G below the 
position of the zero-energy resonance. 

In order to gain further insight into the range of applica- 
bility of the experimental pairwise projection technique, we 
have analysed the critical temperature for the BCS transi- 
tion (depicted in Figure|8} as well as the fermionic condensate 
fraction /c, i.e. the density of condensed pairs divided by one 
half of the total atomic density. This quantity /c is thus given 
by 



/c 



J dp \Kp 



nil 



(56) 



One of the measurable quantities in the pairwise projection 
technique is the density of condensed molecules observed af- 
ter the fast sweep of the magnetic field strength below the 
position of the zero-energy resonance. For an idealised, i.e. 
instantaneous, sweep, this condensate density would be de- 
termined just by the projection of the initial pair wave func- 
tion onto the wave function of the highest excited vibrational 
bound state 0b(5proj) [i-e. an exact molecular bound state of 
the two-channel two-body Hamiltonian ^] at the final mag- 
netic field strength Bpioj. This Feshbach molecular state be- 



FIG. 8: (Color online) The critical temperature for the BCS tran- 
sition as a function of the magnetic field strength B for a range of 
densities n. The density of n = 1.5 x 10''' cm"^^ corresponds to the 
Fermi temperature of Tf = 0.35 yuK quoted in Ref. |9]. We note the 
slight differences between the predictions of the boson-fermion (thin 
solid curve) and standard fermionic (dotted curve) for the highest 
density of n = 5 x lO'"* cm"''. This high-density critical temperature 
curve reflects the only calculation for which we were able to identify 
any visible differences between the two approaches. 



comes populated by the sweep. In the case of Ref. j9|] the final 
magnetic field strength was Bpioj = Bq - 10 G. We have nor- 
malised the density of condensed molecules by a half of the 
atomic density «, in analogy to the definition of the fermionic 
condensate fraction of Eq. ( I56> . This yields: 



/mol 



n/2 



(57) 



where (p^^ is the entrance-channel component of the two- 
channel Feshbach molecular bound state given by Eq. (I16t 
and Nb is the bound-state normalisation constant defined in 
Eq. (I18> . The fermionic condensate fraction /c and its overlap 
/mol with the bound state wave function <^b(fiproj) are shown in 
Fig.|9l(see also (^). 

As shown in Fig. |5] at the lowest experimentally reported 
t9J temperature of 0.08 Tp the fermionic condensate should 
be present up to about 0.6 G above the zero-energy resonance. 
For magnetic field strengths larger than about Bq +0.5 G, how- 
ever, the overlap of the fermionic pair wave function with the 
target molecular wave function is negligible. Consequently, 
the paired fraction at higher magnetic fields is not detected by 
the pairwise projection technique, at least, under the idealis- 
ing assumption that the magnetic field sweep is infinitely fast. 
The magnetic field strength at which the overlap of the paired 
fraction with the target molecular state vanishes also coincides 
with the field strength at which the extent of the pairs equals 
the mean inter-atomic spacing of the gas. 
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FIG. 9: (Color online) The fermionic condensate fraction and its 
overlap /n,oi with the target molecular bound state at 10 G below the 
zero-energy resonance for T = (dashed and solid curves, respec- 
tively) and for T = O.OSrp (dotted and dot-dashed curves, respec- 
tively). The density is « = 1.5 x lO" cm"' corresponding to the 
Fermi temperature of = 0.35yuK quoted in Ref. 



VI. CONCLUSIONS 

We have shown in this paper that the parameters of the 
general standard fermionic and boson-fermion many-body 
Hamiltonians, commonly used in studies of the BCS-BEC 
cross-over problem in dilute atomic gases, can be determined 
in an unambiguous manner by applying the Hamiltonians to 
a single pair of atoms and determining the low-energy bi- 
nary bound state and scattering properties. In particular, we 
have described a convenient two-channel approach to the res- 
onance enhanced two-body scattering, involving separable 
inter-atomic interactions. In a possible extension of this ap- 
proach to the many-body description of a gas the five param- 
eters of the boson-fermion Hamiltonian are all determined 
through five physical parameters of the resonance enhanced 
binary scattering, all of which can be deduced from measure- 
ments of binary scattering properties. 

We have studied in detail the 202 G resonance in a mixture 
of (/ = 9/2, m/ = -9/2) and (/ = 9/2, m/ = -7/2) Zeeman 
states of "'"K fermions and shown that the maximum closed- 
channel admixture in the highest excited vibrational bound 
state reaches only about 8% at 7.2 G below the position of 
the zero-energy resonance. This small admixture is due to an 
avoided crossing between the energies of the Feshbach res- 
onance level and a bound state of the background scattering 
potential. This points to the necessity of including the highest 
excited vibrational state of the background-scattering poten- 
tial, in addition to the background scattering length, in the de- 
scription of binary scattering in the vicinity of the zero-energy 
resonance relevant to the BCS-BEC cross-over Including this 
second parameter characterising the entrance channel scatter- 
ing also significantly extends the range of applicability of ef- 



fective single channel approaches, a fact that often appears 
not to be appreciated. In fact, for the ^*'K system studied in 
our paper we conclude that the energies of the highest-excited 
molecular bound state given by the effective single-channel 
approach are virtually indistinguishable from the results of a 
more elaborate two-channel analysis over the whole range of 
magnetic fields relevant to the experiment of Ref. |9]. The 
universal regime in the vicinity of the zero-energy resonance, 
however, where the binding energy is a function of the scat- 
tering length only, covers only a rather small fraction of the 
BCS-BEC regime in the case of ^^K. 

These observations suggest that the effective single-channel 
approach with only two parameters should be sufficient also 
in studying the many-body BCS-BEC cross-over problem. 
To verify this, we have performed the mean-field analysis 
of the cross-over regime using the standard fermionic and 
the commonly used boson-fermion Hamiltonians for a wide 
range of atomic densities and magnetic fields. As anticipated, 
the results of our comparative studies remain virtually in- 
distinguishable for all the many-body observables common 
to both types of Hamiltonians. The population associated 
with the bosonic field of the boson-fermion Hamiltonian re- 
mains, for the whole range of experimentally relevant densi- 
ties, very close to the closed channel populations of the Fesh- 
bach molecule predicted by the two-body two-channel ap- 
proach. In particular, in the regime where it reaches its maxi- 
mum value of only about 8% (at 7.2 G below Bq) this quantity 
is already virtually density independent. For B > Bq and the 
density of the experiment of Ref. |9] it remains below 0.6%, 
in stark contrast to the results of Refs. L32. .33J . We there- 
fore conclude that for the '^''K system the standard fermionic 
many-body approach should be sufficient to study the whole 
BCS-BEC cross-over regime, as long as it properly accounts 
for the scattering length and the energy of the highest-excited 
vibrational bound state of the background-scattering potential. 

Adjusting the parameters of the many -body Hamiltonian on 
the basis of the two-body considerations puts the proceeding 
analysis of the BCS-BEC cross-over problem of ""'K on firm 
ground. Adopting a different procedure may lead to unphys- 
ical conclusions. For example, one may reach a conclusion 
that below the zero-energy resonance there is a large popula- 
tion associated with the quantum field describing structureless 
bosons, as it appears in the boson-fermion Hamiltonian. One 
can also mistakenly associate this population with the num- 
ber of the actually produced weakly-bound molecules. Such 
conclusions may result purely from ignoring the dominant in- 
fluence of the background scattering potential and its high- 
est excited vibrational bound state of energy E^i. We have 
demonstrated that the background-scattering potential and its 
highest-excited vibrational bound state are crucial to a proper 
description of both two- and many-body physics outside the 
small universal regime in the vicinity of the zero-energy reso- 
nance. 

We would also like to emphasise that the position Bo of 
the zero-energy resonance, where the weakly bound molec- 
ular state emerges, is, due to strong inter-channel coupling, 
very different from the magnetic field strength Bres where the 
energy of the Feshbach resonance level fires crosses the dis- 
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sociation threshold of the entrance channel. In fact, for the 
"^"K resonance at Bq - 202.1 G, Ey^^iBo) is as large as /ix21.8 
MHz and negative. We note that the analyses of Refs. i32ll33ll 
do not capture this important physical fact as they identify Bres 
with Bo- 

Finally, we have analysed thermodynamic quantities asso- 
ciated with the BCS-BEC cross-over in "'"K. We found that for 
densities of the experiments of Ref. |9] the cross-over point, 
as characterised by a qualitative change in the form of the pair 
wave function from oscillatory to decaying and in the form 
of the single-particle excitation spectrum, takes place about 
0.3 G below the position of the zero-energy resonance. We 
have also studied the pairwise projection technique employed 
in the experiment of Ref. |9] by calculating the overlap of the 
fermionic condensate wave function with the target molecu- 
lar state at about 10 G below the zero-energy resonance. We 
found that when the fermionic pairs start to spatially overlap. 



which in the case studied takes place about 0.5 G above the 
position of the zero-energy resonance, the sensitivity of the 
pairwise projection technique decreases. This result coincides 
with the reported measurements of the condensed fraction of 
fermionic pairs extending up to about 0.5 G above the position 
of the zero-energy resonance. 
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